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Abstract 


Traffic  signals  can  be  synchronized  so  that  a  car,  starting  at 
one  end  of  a  main  artery  and  traveling  at  preassigned  speeds,  can 
go  to  the  other  end  without  stopping  for  a  red  light.  The  portion 
of  a  signal  cycle  for  which  this  is  possible  is  called  the  band* 
width  for  that  direction.  Ordinarily  the  bandwidth  in  each  direc¬ 
tion  is  single,  i.e.,  is  not  split  into  two  or  more  intervals  within 
a  cycle.  For  this  case  we  formulate  the  arterial  problem  as  a 
mixed-integer  linear  program:  Given  (1)  an  arbitrary  number  of 
signals,  (2)  the  fraction  of  the  cycle  that  is  red  at  each  signal, 

(3)  upper  and  lower  limits  on  signal  period,  (4)  upper  and  lower 
limits  on  speed  each  way  between  adjacent  signals,  (5)  limits  on 
change  in  speed,  and  (6)  a  constant  of  proportionality  between  the 
two  ba  idwidths,  find  (1)  a  common  signal  period,  (2)  speeds  between 
signals,  and  (3)  the  relative  phasing  of  the  signals  so  as  to  maxi¬ 
mize  the  sum  of  the  bandwidths.  Several  variants  of  the  problem  are 
also  formulated. 

A  branch  and  bound  algorithm  is  developed  for  solving  *-he  given 
mixed-integer  linear  program  by  solving  a  sequence  of  ordinary  linear 
programs.  A  10  signal  example  is  worked  out. 

The  problem  of  synchronizing  a  network  of  signals  is  also  formu¬ 
lated.  The  resulting  program  consists  of  the  arterial  programs  for 
the  individual  streets  plus  a  set  of  further  constraints  that  arise 
because  the  streets  connect  together  to  form  closed  loops.  The  objec¬ 
tive  function  used  for  the  network  program  is  a  weighted  sum  of  the 
bandwidths  in  each  direction  on  each  artery.  A  7  signal  example  is 
worked  out. 
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1.  Introduction 


Traffic  signals  can  be  synchronized  so  that  a  car  can  go  from 
one  end  of  a  street  to  the  other  without  stopping,  if  the  driver 
maintains  the  speeds  used  in  setting  the  signals.  The  portion 
of  the  signal  cycle  for  which  this  is  possible  is  called  the  band¬ 
width  for  that  direction.  Synchronizations  far  large  bandwidth  are 
frequently  called  progressions.  Traffic  engineers  have  long  set  up 
progressions  and  they  seem  to  be  quite  helpful  when  traffic  is 
light,  although  possibly  of  little  value  when  traffic  is  heavy. 

The  work  of  Helly  and  Baker  ^  tends  to  suggest  that  this  is  the 

case.  The  value  of  a  progression  can  be  enhanced  by  giving  the 

[2] 

drivers  speed  instructions,  as  has  been  done  by  von  Stein. 

[3] 

In  an  earlier  paper  J.  Morgan  and  the  author  have  developed 
an  algorithm  for  solving  the  following  two  problems  of  synchronizing 
traffic  signals  for  a  progression  on  an  arterial  street,  in  the 
case  where  bandwidths  in  each  direction  are  unsplit. 

1.  Given  an  arbitrary  nuotber  of  signals,  a  common  signal 
period,  the  green  and  red  times  for  each  signal,  and  speci¬ 
fied  travel  times  between  adjacent  signals,  synchronise 
the  signals  to  produce  bandwidths  that  are  equal  in  each 
direction  and  as  large  as  possible. 

2,  Adjust  the  synchronisation  to  Increase  one  bandwidth  to 
some  specified,  feasible  value  and  maintain  the  other  as 
large  as  is  then  possible. 

Subsequently,  S.  Oliver  pointed  out  to  the  author  that  a  two-signal 
version  of  Problem  1  could  be  set  up  as  a  linear  program.  We  here 
expend  this  idee  into  a  rather  general  formulation  of  the  maximal 
bandwidth  problem.  In  the  general  case  we  still  have  a  linear 
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program  but  unfortunately  it  is  of  the  mixed-integer  type.  For 
Problems  1  and  2  the  new  formulation  offers  no  advantages  and  many 
disadvantages.  However,  the  linear  programming  format  opens  up  the 
possibility  of  solving  more  general  problems  involving  the  intro¬ 
duction  of  new  decision  variables  and  new  constraints. 

For  example,  maximal  bandwidth  calculations  as  usually  performed 
have  a  disconcerting  feature.  On  a  long  street  the  critical  signals 
that  constrict  bandwidth  may  turn  out  to  be  very  far  apart.  Then 
a  small  change  in  one  of  the  design  speeds  along  the  street  is 
likely,  upon  re-solving  the  problem,  to  result  in  a  different  syn¬ 
chronization  and  a  different  bandwidth.  Yet,  drivers  do  not  hojd 
their  speeds  exactly  constant  and,  indeed,  as  is  well  known,  tend 
to  adjust  their  speeds  to  the  signals.  Thus,  it  makes  a  good  deal 
of  sense  to  let  the  design  speed  between  signals  be  a  variable,  at 
least,  within  certain  limits.  This  can  easily  be  done  in  the  linear 
program. 

Another  variable  that  can  be  introduced  explicitly  is  the 
signal  period.  In  Problem  1,  period  is  a  constant  and,  although 
it  is  not  too  difficult  to  use  our  earlier  methods  to  examine  a 
considerable  number  of  values  in  an  organised  way,  the  continuous 
variation  of  the  linear  program  formulation  seams  preferable. 

Perhaps  the  most  Interesting  development  is  that  the  problem 
of  synchronising  signals  for  a  network  of  streets  can  be  formulated 
as  a  mixed-integer  linear  program.  The  program  for  a  network  con¬ 
sists  of  arterial  programs  for  individual  streets  plus  additional 
constraints  whenever  the  arteries  connect  together  to  form  loops  or 
cycles. 
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2.  The  Arterial  Problem 

2.1  Definitions .  Consider  a  two-way  street  having  n  traffic 
signals.  Directions  on  the  street  will  be  identified  as  outbound 
and  inbound.  The  signals  will  be  denoted  S^,  S  ,...,Sn  with  the 
subscript  increasing  in  the  outbound  direction. 

Figure  1  shows  a  space-time  diagram  for  travel  on  a  street. 

Heavy  horizontal  lines  indicate  when  the  signals  are  red.  The  zig¬ 
zag  lines  represent  tra lectories  of  the  cars  passing  unimpeded  along 
the  street  in  the  directions  shown.  Changes  in  slope  correspond  to 
changes  in  speed.  The  set  of  possible  unimpeded  trajectories  in  a 
given  direction  forms  a  green  band  whose  horizontal  width  is  the 
bandwidth  for  that  direction.  Although  drawn  but  once,  the  green 
bands  occur  once  per  cycle  in  parallel  bands  across  the  diagram. 

Ordinarily  the  bandwidth  in  each  cycle  is  single,  i.e.,  is  not 
split  into  two  or  more  intervals  within  a  cycle.  A  split  bandwidth 
can,  however,  occur  and  special  examples  have  been  constructed  in 

which  the  maximal  total  bandwidth  is  made  up  of  two  pieces  in  one  or  both 

[A] 

directions.  This  possibility  will  be  ignored  in  the  present  formu¬ 
lation..  The  mathematical  program  to  be  constructed  will  maximize 
the  weighted  sum  of  twu  bandwidths,  one.  taken  from  each  direction 
without  considering  the  possibility  that  other  pieces  may  exist. 

Curtain  signals  form  the  uiti  ite  limitation  of  bandwidth  and 
will  be  called  critical  signals.  A  signal  Sj  is  said  to  be  a 
critical  signal  if  one  side  of  Sj's  rad  touches  the  green  band  in 
one  direction  and  the  other  side  touches  the  green  band  in  the 
other  direction.  Thus,  in  Figure  1  signals  and  are  critical, 
but  no  others  are. 
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Time 

(cycles) 


2.2  Basic  Maximal  Bandwidth  Formulation.  First  we  set  up 
Problem  1  above  as  a  mixed-integer  linear  problem.  Let 

r^  =  red  time  of  on  street  under  study,  (cycles) 

b  (b)  =  outbound  (inbound)  bandwidth,  (cycles) 

t(h,i)  (t(h,i))  =  travel  time  from  to  in  outbound  direc¬ 
tion  (travel  time  from  to  in  inbound 

direction).  These  quantities  are  algebraically 
positive  if  S.  follows  in  the  outbound 
direction,  otherwise  negative,  (cycles) 

0(h,i)  (0(h,i))  =  distance  from  center  of  red  at  to  the 

center  of  a  particular  red  at  S^.  See  Figure 
2.  The  two  reds  are  chosen  so  that  each  is 
immediately  to  the  left  (right)  of  the  same 
outbound  '(inbound)  green  band.  0(h,  i)  (/(h,i)) 
is  positive  if  S^s  center  of  red  lies  to  the 
right  (left)  of  Sh's.  (cycles) 

Wi  ^Wi^  =  distance  ^rom  the  rifcht  (left)  side  of  S^'s 

red  to  the  green  band.  See  Figure  2.  (cycles) 
m(h,i)  =  0(h,i)  +  0(h,i)  (2.1) 

Note  that  a  quantity  having  the  dimensions  of  time  can  always  be 
expressed  in  cycles  by  dividing  by  the  period. 

From  Figure  2  can  be  read  the  identities: 

(1/2) ru  +  wL  +  t(h,i)  -  w.  -  (1/2) r.  =  0(h,i)  (2.2a) 

h  n  i  * 

(l/?.)rh  + w^+  t(h,i)  -  w^  -  (i/2)ri  =  0(h,i)  .  (2.2b) 
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The  values  of  0(h, i)  and  0(h, i)  have  the  very  important  constraint 
that  their  sum  must  be  an  integer,  but  otherwise  are  unrestricted. 
Therefore,  the  above  two  constraints  can  be  replaced  by 


(wh+  wh)  -  (wt+  +  t(h,i)  +  t(h, i) 


m(h, i)  =  integer  . 


m(h,  i)  -  (rh-  r^ 
(2.3) 


From  physical  considerations,  it  follows  that 


t(h,j)  =  t(h,i)  +  t(i,j) 
t(h,i)  =  -  t(i,h) 

Therefore, 

0(h,j)  =  0(h,i)  +  0<i,j) 

0(h,k)  =  -  0(i,h) 


(2.4a) 


(2.5a) 


(2.4b) 

(2.5b) 


The  corresponding  inbound  expressions  are  obtained  by  adding  bars 
to  (2.4)  and  (2.5)  above.  Then, 

m(h,J)  =  m(h,i)  +  ra(i,j)  (2.4c) 

m(h,i)  =  -  m(i,h)  .  (2.5c) 

From  Figure  2  we  also  see  that 

*1  +  b  <  l  -  r  (2.6a) 

w1  +  b  <  1  -  r£  .  (2.6b) 
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So  far,  the  expressions  apply  to  any  arbitrarily  indexed 
signals.  No  use  has  been  made  of  the  special  numbering  along  the 
outbound  direction.  In  later  work  on  networks  the  present  gener¬ 
ality  will  be  needed,  but  we  can  simplify  notation  here  by  de¬ 
fining 

xt  =  x(i,  i  +  1)  (2.7) 

for  x  =  t,  t,  0,  0,  m  . 

Problem  1  will  now  be  represented  by  a  mixed-integer  linear 
program.  Constraints  (2.6)  are  needed  for  each  signal  and  (2.3) 
for  each  adjacent  pair.  (Although  there  are  many  other  types  of 
pairs,  it  suffices  to  consider  adjacent  pairs  since  the  constraint 
(2.3)  for  any  other  pair  can  be  obtained  from  linear  combinations 
of  (2.3)  for  adjacent  pairs.) 

LP1 .  Find  b,  b,  w w ,  m  to 

v/vv  XXX 

max  (b  +  b)  , 


Subject  to: 


b  =  b 

(LPl.i) 

+  b  <  1  ’  ri 

(LP1.2e) 

y  i»l,...,n 

+  F  <  1  -  1 

(LP1.2b) 

<wt+  v  -  (wi+i+  Vi> +  (ti+  V  =  Bi  ■  (v  ri+i)  <lpi-3> 


1*1 , . • • , n* 1 


integer 
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t>,  b>  vt>  wt  >  0 


LP1  has  3n  equations  and  3rvf  1  unknowns,  not  counting  slack  or 
artificial  variables. 

2.3  Adding  Decision  Variables.  Next  we  let  period  and  speed 
be  variables.  Each  will  be  constrained  by  upper  and  lower  limits. 
In  addition,  changes  in  speed  from  one  street  segment  to  the  next 
will  be  limited.  Finally,  instead  of  requiring  equal  bandwidths 
in  each  direction,  we  let  one  bandwidth  be  a  fixed  proportion  of 
the  other.  (However,  this  is  Just  one  of  the  possible  constrain¬ 
ing  relations  that  might  be  used.) 

Let 

k  =  constant  of  proportionality  between  b  and  b. 

T  =  signal  period,  (seconds) 

Tj^,  ~  l°wer  and  upper  limits  on  period,  i.e., 

<  T  £  T2  *  (seconds) 
z  a  signal  frequency,  (cycles/second) 

d(h,i)  a  distance  between  and  S^.  (meters) 

dt  a  d(i,  i+1) 

vi  ^i}  "  8P*«d  between  and  outbound  (Si+1  and 

inbound) .  (meters/second) 

ei'^i  ^i'^i^  *  lover  and  upper  limits  on  outbound  (inbound) 
speed,  (meters/second) 

l/h^l/g^Cl/K^jl/gj,)  =  lover  and  upper  limits  on  change  in  outbound 

(inbound)  reciprocal  speed,  i.e. 
l/ht  <  U/v1+l)  -  d/vt)  <  l/g£  . 

-1 


(meters/second) 
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We  are  constraining  change  in  speed  by  putting  upper  and  lower 
limits  on  change  in  reciprocal  speed.  Although  the  two  are  not 
quite  the  same,  constraining  the  change  in  reciprocal  speed  surely 
satisfies  the  basic  intention,  which  is  to  have  a  means  of  pre¬ 
venting  large  abrupt  speed  changes.  The  reason  for  working  with 
reciprocal  speed  is  that  it  appears  linearly  in  the  constraints  and 
can  be  transformed  into  t  .  Thus 

tf  =  (d./v.)z  ,  t.  =  (di/v.)z  .  (2.8) 

In  an  expanded  program  t^,t^,  and  z  are  decision  variables,  which, 
once  known,  determine  the  speeds. 

After  some  algebraic  manipulations  the  constrainsts  on  band¬ 
width,  period,  speed,  and  change  in  reciprocal  speed  can  be  added 
to  LPl  to  give  a  more  versatile  mixed-integer  linear  program. 

LP^.  Find  b,  b,  z,  w^.,  t^  t^  to 
max  (b  +  b) 

subject  to: 


b  =  k  b 

(LP2.1) 

1/T2  <  : 

!  <  1/T1 

(LP2.2) 

+  b  < 

*•'*  1 

c 

*% 

• 

• 

• 

II 

•H 

(LP2.3a) 

wi  +  F  < 

l-*t  J 

(LP2.3b) 

(V-V-(”i+i*W+(ti+Ti>  *  v(trri*ih  <LM-4> 

m^  a  integer 

(di/fi)z 

*  *1  S 

/  (LP2.Se) 

VI 

|U 

VI 

(di/*i)* 

|  (LP2.5b) 

-10- 


(dt/h^z  <  (dt/di+1)ti+i  -  tt  <  (di/g.)z  j  (LP2.6a) 

}  n-2 

(di/Vz  ^  (di/di+i)ti+i  •  Ei  ^  (di/ii)Z  j 


(LP2. 6b) 


J 


b,  b,  w.  >  0 
1  ’  i  l  — 


LP2  involves  (lln-10)  constraints  and  5n  variables,  not 
counting  slack  and  artificial  variables. 

2.4  Symmetric  Problems.  Certain  symmetries,  if  they  exist, 
can  be  used  to  reduce  the  size  of  LP2.  In  particular,  if  the  con¬ 
straints  on  speed  are  the  same  in  each  direction,  and  if  bandwidths 
in  each  direction  are  required  to  be  equal,  the  size  of  the  linear 
program  is  cut  about  407.. 

A  time-saving  technique  may  be  to  solve  for  maximal  equal 
bandwidths  and  then  adjust  the  synchronization  for  other  cases. 
Reference  3  shows  that,  with  some  qualifications  to  cover  special 
conditions,  max  (b  +  b)  is  a  constant  that  can  be  divided  between 
the  directions.  Rules  are  given  there  for  setting  one  bandwidth 
to  any  feasible  value  and  the  other  to  the  largest  value  that  is 
then  possible  for  the  given  speeds  and  period. 

The  linear  program  for  the  symmetric  case  is  based  on  the 
following  two  theorems. 

Theorem  1.  If  the  limits  on  speed  and  change  in  speed  are 
the  same  in  each  direction,  and  If  LP2  has  any  optimal  solutions, 
then  LP2  has  an  optimal  solution  in  which  t^  =  t^  for  each  i. 

Proof.  To  hove  the  same  limits  on  speed  and  change  in  speed 
in  each  direction  is  to  have  eA  <=  e^,  =  f^,  gi  =  g^,  and  h^  =  h^. 
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Let  j t ti|  be  the  set  of  of  t's  in  an  optimal  solution  for  LP2. 

Define 

t  *  =  t  1  =  (ti  +  tt)/2  ,  i  =  l,...,n-l  . 

We  claim  that  the  t^  t^1  and  t^  <r  t^'  will  yield  a  new  optimal 
solution  for  LP2.  Since  the  change  of  t's  does  not  affect  the  ob¬ 
jective  function,  the  only  question  is  feasibility.  Clearly, 

LP2.4  is  still  satisfied.  Adding  LP2.5a  to  LP2.5b  and  dividing  by 
2  gives 

(di/fi)z  <  tt'  <  (di/ei)a  , 

so  that  LP2.5a  and  b  are  satisfied.  A  similar  argument  shows  that 
LP2.6a  and  b  are  satisfied.  This  completes  the  proof. 

Theorem  2.  If  the  bandwidths  in  each  direction  are  required 
to  be  equal,  and  if  LP2  has  any  optimal  solutions,  then  LP2  has  an 
optimal  solution  in  which  w^  =  for  each  i. 

Proof.  Given  any  optimal  solution  to  LP2  with  b  =  b  ,  replace 
and  w^  by  w^'  and  w  '  as  defined  by 

wi'  •  wi’  =  (wt  +  wi)/2  . 

The  same  type  of  argumenti  used  to  Theorem  1  demonstrate  the  fees- 

t 

ibility  of  the  new  solution  in  LP2.3a  and  b  and  LP2.4.  Optimality 
is  then  immediate  since  the  objective  function  is  unchanged. 

If  the  conditions  of  both  theorems  are  satisfied  we  can  formu¬ 
late  the  following  simplified  mixed-integer  linear  programs 

LP3.  Find  b,  a,  w. ,  t.,  m.  to 
111 

max  b 
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*9 


Subject  to; 


(1/T2)  <  z  <  (1/T^ 


w^  +  b  <  1  -  r^ 


i=l, . .  .,n 


VV  +  ci  -  (1/2)”i  -  <l/»<Vri+i> 


mi  =  integer 


i=l, . . .,n-l 


(LP3.1) 

(LP3.2) 

(LP3.3) 


(di/fi>z  -  fci  ~  <Vei)z  ' 

(di/hi)z  ^  (di/di+l)t:i+l  "  Ci  -  (di/8i)Z|  ^  ^ 
b,  >  0 


(LP3.4) 

(LP3.5) 


LP3  has  (6n-S)  constraints  and  3n  variables  not  counting  slack 
and  artificial  variables. 

2.5  Determining  the  Synchronization.  The  linear  programming 
variables  determine  the  synchronisation  of  the  signals.  Let 

6(h,i)  =  relative  phase  (offset)  of  S^  and  S^,  measured 

as  the  time  from  the  center  of  a  red  of  S.  to 

n 

the  next  center  of  red  of  S^.  (cycles) 

An  illustration  of  0<h,i)  is  given  in  Figure  2.  We  adopt  the  con* 
vent ion  0  <  0(h,i)  <  1.  A  set  of  0(1,1),  1=2,..., n  will  be 

called  a  synchronisation  of  the  signals. 

In  order  to  give  simple,  explicit  expressions  for  0  and 
certain  other  quantities,  we  define,  for  arbitrary  real  x: 
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int[x] 
man[ x] 


largest  integer  <  x 
x  -  int[xj 


(2.9) 


(2.10) 

Thus,  for  example,  int  [5.2]  =  5,  man  [5.2]  =  .2,  int  [-1.3]  =  -2, 
man  [-1.3]  =  .7  . 

From  Figure  2  and  (2.2)  we  see  that 

9(h, i)  =  man[0(h,i)]  (2.11) 

i-1 

6(1, i)  =  man[w.-w  +  £  fck+(l/2)  (r, -r  )  ]  (2.12) 

L  1  k=l 

In  the  symmetric  case  of  LP3,  (2.2)  shows  that  0(h,i)  =  0(h, i), 
whence 


0(h,i)  =  (l/2)m(h, i) 

(2.13) 

i-1 

, 

6(1,  i)  =  man[  (1/2)  £  m.]  . 

k=i  K 

(2.14) 

2.6  Limits  on  the  integer  Variables. 

The  integer  variable 

in(h,i)  is  rather  limited  in  the  values  it  can  take  on  without 
causing  an  infeasible  program.  If  most  of  the  infeasible  values 

can  be  excluded  in  advance,  wr.  should  be  able  to  save  computation 

% 

time. in  the  algorithm  of  the  next  section. 

He  seek  integers m'(h,i)  and  m"(h,i)  such  that 
m'(h,i)  <  m(h,  i)  <  mf\h,i)  . 

From  (2.3) 

m(h,i)  =.  <wh+wh)  ”  (wi*^l>  +  t<M)  +  t(h,i)  +  (rh‘ri^  * 

Since  0  <  <  1  -  r^  for  any  k  and  since  a  similar  relation 

holds  for  w^  , 
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ft. 


t(h,i)+  t(h,i)-(l-rh)-(l-ri)  <  m(h,i)  <  t(h,i)+  t(h,i)+  (l-rh)+(l-rt)  . 


Therefore, 

nKh,i)  =  int[(l-rh)  +  (1-r.)  +  t(h,i),+  t(h,i)J  (2.15a) 

m'(h,i)  =  -int[(l-rh)  +  (1-r^  -  t(h,i)  -  t(h,i)]  (2.15b) 


However,  t  and  t  are  usually  variables  and  can  be  bounded  in  terms 
of  the  constants  of  the  problem.  Thus,  in  LP2, 


VfiT2  -  VS  Wi 


Following  the  notation  of  LP2,  we  obtain: 


m"  . 
i 


intL(l-rt)  +  (l-ri+1)  +  (d^e^)  +  (d^e^)]  (2.16a) 


=  -intKl-r^  +  (1-r^)  -  (d^f^)  -  (\fiff].  (2.16b) 

3.  A  Branch  and  Bound  Algorithm 

The  algorithms  presently  available  for  solving  the  general  mixed- 
integer  linear  programming  problem  are  usually  considered  only 
partially  satisfactory.  Relatively  small  problems  sometimes  take 
an  inordinate  length  of  time  to  solve.  Consequently,  a  specialized 
algorithm  has  been  developed  for  the  maximal  bandwidth  problem  by 
using  branch  and  bound  methods.  These  methods  have  proven  rather 
successful  in  some  other  combinatorial  problems. 

3.1  General  Description.  The  basic  idea  of  the  method  is  to 
break  up  the  set  of  all  feasible  solutions  into  smaller  and  smaller 
subsets  and  to  calculate  for  each  an  upper  (lower)  bound  on  the 
objective  function  of  the  best  solution  therein.  The  bounds  guide 
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the  partitioning  of  the  subsets  and  eventually  identify  a  maximizing 
(minimizing)  solution.  When  a  subset  is  found  that  contains  a  single 
solution  whose  objective  function  is  greater  (less)  than  or  equal  to 
the  upper  (lower)  bounds  for  all  other  subsets,  that  solution  is 
optimal.  The  subsets  of  solutions  are  conveniently  represented  as 
the  nodes  on  a  tree  and  the  process  of  partitioning  as  a  branching 
of  the  tree;  hence  the  name,  "branch  and  bound". 

The  convergence  of  the  process  can  be  assured  by  devising  a 
partitioning  procedure  that,  at  worst,  leads  to  eventual  enumeration 
of  all  solutions.  The  computational  efficiency  of  the  process,  how¬ 
ever,  is  very  dependent  on  the  methods  used  to  perform  the  partitions 
and  calculate  the  bounds.  The  successful  applications  to  date  have 
exploited  special  features  of  the  particular  class  of  problems  at 
hand  to  develop  good  bounds. 

A  branch  and  bound  algorithm  for  LP2  (or  LP3)»/ill  now  be  de¬ 
scribed.  Consider  an  ^-signal  problem  (2  <  r  <  n)  consisting  of 
the  first  r  signals  of  the  n-signal  problem.  Suppose  that  we  pick 
specific  integer  values  for  the  variables  m^,...,mr_^  *n  LP2.  LP2  is 
now  an  ordinary  linear  program  and  can  be  solved  in  a  straight¬ 
forward  manner  to  maximise  the  objective  function.  Now  the  set  of 
n-signal  solutions  that  have  the  specific  values  m^...mr_^  is  a 
subset  of  the  set  of  all  n-signal  solutions.  Furthermore,  the 
maximal  objective  function  for  the  r-signal  problem  having  these 
m-values  forms  an  upper  bound  on  the  objective  function  for  any 
n-signal  solution  with  these  m-values,  since  adding  more  signals  can 
only  restrict  the  possibilities  for  maximisation.  Thus  we  have  a 
method  for  defining  subsets  of  solutions  (by  specifying  m-values) 
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and  for  placing  bounds  on  their  objective  functions  (by  solving 
ordinary  linear  programs) .  The  subsets  can  be  further  partitioned 
by  specifying  values  for  more  of  the  m's. 

3.2  Flow  Diagram.  The  ideas  of  the  previous  section  will  now 
be  expressed  more  formally.  Let 

X  =  a  subset  of  solutions  of  the  n-signal  problem; 

namely,  the  set  of  all  solutions  that  have  the  m- 
values  of  a  particular  r-signal  problem. 
r(X)  =  the  number  of  signals  in  the  problem  defining  X. 
m^X), . .  .,m  ^(X)  =  the  integer  values  of  m^  for  the  r-signal 

problem  defining  X. 

Xj ,  X^...  =  subsets  of  solutions  obtained  by  partitioning  X 

further. 

The  set  X  is  partitioned  into  subsets  x^>  X0, ...  each  defined  by 
a  different  (r+1) -signal  problem.  Any  of  these  (r+1) -signal  prob¬ 
lems,  say  the  one  for  X^,  will  have  the  same  m^'s  as  X  up  to 
i  =  r-1: 

=  mi(X)  i=l,...r-l  , 

and,  in  addition,  a  value  for  m^.  Let 

m  m  ”  =  smallest  and  largest  value  of  m  to  be  used  in 

r  r  w  r 

partitioning  X. 

a(X)  =  an  upper  bound  on  the  objective  function  for  the 

solutions  in  X,  as  obtained  by  solving  an  ordinary 
linear  program  for  the  r-signal  problem  with 
specified  m-values. 
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0  =  best  n-signal  objective  function  found  so  far  in 

the  calculation. 

A  flow  diagram  of  the  algorithm  is  shown  in  Figure  3.  We 
shall  trace  through  the  boxes.  As  this  is  done,  it  may  be  helpful 
to  refer  to  1  jure  4,  which  shows  the  tree  of  a  worked-out  10 
signal  problem.  In  the  tree  representation,  subsets  of  solutions 
are  represented  by  nodes  and  so  we  shall  use  X  to  refer  to  both 
a  subset  and  its  node. 

Box  1  iuitialires  three  quantities:  X  to  the  set  of  all 
solutions,  r  to  l,  and  0  to  zero.  Box  2  initializes  two  quanti¬ 
ties  to  be  used  in  the  construction  of  the  nodes  branching  out  of 
X:  The  index  .j  for  counting  the  nodes  is  set  to  1  and  m^,  whose 

values  will  define  the  nodes,  is  set  to  its  lower  limit  m  '. 

’  r 

Box  3  sets  up  a  node  Xj  on  a  branch  out  of  X.  The  node  may 

be  identified  by  m  (X.),  which  ia  set  to  the  current  value  of  m  . 

r  j  r 

(In  Figure  4  the  value  of  m^(X^)  is .written  inside  the  circle  of 
each  node,  X^>  To  define  the  node  completely,  we  also  need 
m1(Xj),...,mr_1(Xj),  but  these  can  be  found  by  tracing  back 
through  the  tree  to  "all  solutions"  and  reading  off  the  m^  in  the 
nodes  encountered.  Next  LP2  is  solved  for  the  (pfl) -signal  prob¬ 
lem  defined  by  the  current  m-values,  m^, ...,m^.  The  maximal 
value  of  the  objective  function  becomes  the  upper  bound  a(X^)  and 
Is  used  to  label  the  node. 

Box  4  tests  whether  an  n-signal  problem  is  currently  being 
solved.  If  so,  box  5  tests  to  see  whether  the  resulting  objec¬ 
tive  function  is  better  than  the  previous  best.  If  this  is  the 
case,  box  6  replaces  the  previous  best  with  the  new  solution. 
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Search  the  terminal  nodes 
having  r  <  n  for  a  node 
with  the  largest  value 
of  a.  Set  X  to  be  this 
node.  Set: 

r * r (X) 

mj«m,  (X)  i»l,..,  r-l 


10  a  (X)  >  £  ? 

'  in° 

Finish 

Last  saved  solution 
is  optimal. 


Figure  3.  Flow  diagram  of  a  branch  and 
bound  algorithm  for  LF2. 
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Box  7  tests  to  see  whether  all  the  planned  mr  values  have 

been  used  to  create  new  nodes.  If  not,  box  8  increases  j  and 

by  one  and  returns  to  box  3  to  set  up  a  new  node. 

Box  9  searches  for  a  new  node  from  which  to  branch.  The 

candidates  are  the  terminal  nodes  of  the  tree,  i.e.,  those  with¬ 
out  branches.  There  is  no  need  to  consider  nodes  with  r  =  n. 

Of  the  eligible  nodes,  the  one  with  the  largest  a  is  selected  and 
becomes  the  new  X.  The  current  value  of  r  is  reset  to  r(X)  and 
the  current  values  of  to  m^(X)  for  i=l,...,r-l. 

Box  10  tests  to  see  whether  the  largest  upper  bound,  a,  is 
greater  than  the  best  n-signal  objective  function  so  far.  If  so, 
there  is  some  possibility  of  finding  a  larger  n-signal  objective 
function  and  the  calculation  returns  to  hox  2  to  branch  from  the 
new  X.  Otherwise,  the  calculation  is  finished,  and  the  current 
best  solution  is  optimal  for  the  original  mixed-integer  n-signal 
problem. 

It  may  be  helpful  to  trace  through  the  first  few  nodes  in 
Figure  4.  First  the  node  "all  solutions"  is  laid  out.  Then, 
since  m^'  =  0  and  m^"  =  1,  the  nodes  m^  =  0  with  a  =  .440  and 
=  1  with  a  =  .292  are  laid  out.  The  largest  a  of  these  two 
is  .440  and  so  the  ra^  *  0  node  is  for  branching.  Out  of  it  comes 
the  nodes:  m^  =  0,  a  =  .  282  and  n>2  *  1,  a  *>  .  381.  The  next  step 
would  be  to  branch  from  the  m2  «  l,  a  =  .381  node.  Eventually 
the  tree  is  developed  as  shown.  Notice  that  the  10-signal  node 
labeled  "optimum"  has  an  objective  function  of  .282,  which  is 
greater  than  or  equal  to  the  upper  bound  for  any  of  the  other 
terminal  nodes. 
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.440 


SOLUTIONS 


,282\  .020  .303/3091  \.  1 9 1 


.158  /  .292' 


1531 .261  \.293\  .094 


.293  .214, 


.280 \  .I4t 


x/.232\  \.206  X/284\ 


X/  .261 


.266 


0  I 


-I)  ^0 


.282  .158 


.057 


optimum 


Figure  4.  Tree  for  10  signal 
arterial  problem.  Labels  on 
nodes  are  the  calculated  upper 
bounds  on  bandwidth  (in  cycles) 
for  the  solutions  represented 
by  the  node.  The  label  X  in¬ 
dicates  no  feasible  solution. 
The  numbers  Inside  the  nodes 
are  values  for  the  mi  shown  at 
the  right. 


3.3  Example.  The  10  signal  problem  being  used  as  a  numeri¬ 
cal  example  is  taken  from  Reference  3  and  represents  a  stretch  of 
Euclid  Avenue  in  Cleveland.  The  distances  and  green  splits  are 
taken  from  actual  usage.  The  limits  on  period,  speeds,  and  speed 
changes  have  been  made  up  fiirly  arbitrarily. 

The  constants  of  the  problem  are  briefly  as  follows.  The  sign 
nals,  starting  with  and  working  outbound,  are  located  at  0,  168, 
381,  716,  929,  1173,  1371,  1493,  1706,  and  1843  meters.  The  cor¬ 
responding  red  times  are  .47,  .40,  .40,  .47,  .48,  .42,  .40,  .40,  .40, 
and  .42  cycles.  Limits  on  period  are  T^  =  55  seconds  and  T^  =  75 
seconds.  Lower  and  upper  limits  on  speed  are  e^  =  e^  =  13.4  meters/ 
second  (30  mph,  48.3  km/hr)  and  f^  =  f^  =  17.9  meters /second  (40  mph, 
64.4  km/hr)  for  each  i.  Limits  on  change  in  reciprocal  speed  are 
1/g^  =  1/g^  =  -  1/h^  =  s  1/h^  =  .0121  (meters/second)  1  for  each  i. 
This  corresponds  to  a  maximum  possible  change  in  speed  of  about  ±2.2 
meters/second  (4.9  mph,  7.9  km/hr)  at  the  lower  limit  of  speed  and 
±3.9  meters/second  (8.7  mph,  14.0  km/hr)  at  the  upper  limit  of  speed. 

The  problem  can  be  solved  by  LP2  or, .since  we  have  chosen  to 
solve  for  equal  bandwldths,  by  LP3.  The  tree  developed  by  the  branch 
and  bound  algorithm  appears  in  Figure  4.  A  total  of  52  linear  pro¬ 
grams  are  required.  The  maximal  bandwidth  for  each  direction  is  .282 
cycles.  The  space-time  diagram  for  the  street  is  shown  in  Figure  5. 
The  optimal  period  is  75  seconds,  making  the  bandwidth  21.2  seconds. 
Starting  between  and  S2  « nd  working  up,  the  ppeeds  are:  17.9, 
17.9,  17.1,  14.2,  13.4,  14.9,  13.4,  15.6,  17.9  meters/second.  The 
phases,  8(1,1),  in  order  of  increasing  index  starting  with  i«l, 
are:  0,  0,  0,  1/2,  1/2,  0,  0,  0,  1/2,  1/2  cycles.  The  critical 
signals  are  8^  S^,  S5,  S&,  and  S9  . 
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:e  S.  Space-time  diagram  for  symmetric  10  signal  problem.  Maximal  bandwidth 
seconds.  The  signal  period  is  75  seconds. 


The  results  can  be  compared  to  those  obtained  by  exploring 
various  constant  speeds  and  periods  with  the  methods  of  Reference 
3.  The  best  result  found  was  15.2  meters/second  at  65  seconds 
for  a  bandwidth  of  .235  cycles.  This  is  essentially  the  syn¬ 
chronization  displayed  in  Reference  3.  The  additional  flexibility 
of  variable  speed  has  permitted  at  20%  increase  in  bandwidth. 

3.4  Discussion.  For  the  numerical  example  reported  here,  the 
steps  of  the  branch  and  bound  algorithmwere  carried  out  manually 
except  for  the  linear  programs.  In  some  linear  programming  codes, 
it  is  possible  to  suppress  constraints.  Then  the  entire  n-signal 
problem  can  be  loaded  and  only  those  constraints  corresponding  to 
the  current  r-signal  problem  used.  Changes  in  the  m^  are  ef¬ 
fected  by  changes  in  the  constants  vector.  Frequently,  time  can 
be  saved  by  starting  a  new.  problem  from  the  basis  of  a  similar 
problem  just  completed. 

The  computational  limits  of  the  algorithm  are  relatively  un¬ 
explored.  If  experience  from  a  rather  different  application  ^ 
is  a  guide,  the  number  of  linear  programs  required  is  likely  to 
increase  exponentially  with  n.  If  so,  there  may  well  be  a  fairly 
sharp  upper  limit  to  the  size  of  the  problem  that  is  computation¬ 
ally  feasible.  We  feel,  however,  that  the  example  here  demon¬ 
strates  that  problems  of  a  size  of  practical  interest  can  be  solved. 

A  solution  that  is  likely  to  be  good  but  may  not  be  optimal 
can  be  obtained  by  solving  roughly  2n  linear  programs.  The  pro¬ 
cedure  is  simply  to  branch  next  from  the  node  with  the  largest  a 
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among  those  just  created  from  the  current  X.  Frequently  there 
are  only  two  such  nodes.  The  scanning  of  the  other  terminal 
nodes  is  omitted.  The  value  of  r  increases  monotonically  to  n. 

In  the  example  of  Figure  4,  the  solution  obtained  in  this  manner 
has  9  bandwidth  of  .-278. 

4.  Network  Problems 

The  mixed-integer  formulation  can  be  extended  to  networks. 
The  network  program  consists  of  arterial  programs  for  the  in¬ 
dividual  streets  plus  "cycle  constraints"  that  connect  the 
streets  together.  An  objective  function  can  be  formed  from  the 
bandwidths  of  the  arteries.  A  new  decision  variable,  the  red- 
green  split,  can  usefully  be  introduced  at  certain  signals. 

Some  of  the  decision  variables  of  LP2  will  be  left  out  of 
the  network  formulation  to  be  developed  in  this  section.  The 
reason  is  that  we  can  expect  network  problems  to  be  larger  and 
computationally  more  difficult  than  individual  arterial  problems. 
Therefore  it  seems  desirable  to  show  ways  in  which  the  program 
might  be  simplified.  We  shall  keep  speed  fixed  over  the  length 
of  an  artery,  but  variable  both  for.  the  prtery  as  a  whale  and 
from  artery  to  artery. 

4.1  Arterial  Constraints.  In  a  network  problem  the  term 
artery  will  be  used  to  refer  to  any  street  in  which  a  progression 
is  desired,  i.e.,  on  which  we  define  a  bandwidth  and  bring  it 
into  the  mathematical  program.  As  aarlier,  tha  signals  of  the 
network  are  presumed  to  ha  deaignated  S^,  S2',**'Sn  but  we  can 
no  longar  suppose  that  signals  with  adjacent  index  values  ara 
adjacent.  Consider  than  an  artery  and  supposu  that  tha  signals 
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at  its  extreme  ends  are  and  S^.  We  designate  one  direction, 

say  to  Sy  and  outbound  and  thereby  define  artery  ij. 

Quantities  relating  to  artery  ij  will  be  identified  by  a 
superscript  ij.  Thus,  in  a  straightforward  way,  the  appropriate 
definitions  of  Section  2  generalize  to  b*^,  r^*^,  and  w^^  .  In 
addition,  let 


speed  in  outbound  direction  on  artery  ij. 
(meters/second) 

lower  and  upper  limits  on  outbound  speed  on 
artery  ij.  (meters/second) 

1/v^T  -  reciprocal  speed  outbound  on  artery 

ij.  (meters/cycle) 


Analogous  notation  is  used  as  appropriate  for  inbound  variables. 
When  two  signals  are  required  to  define  a  quantity,  the  arterial 
designation  is  redundant  and  will  be  omitted.  For  example,  the 


quantities  0(h,i), 

0(h,i),  m(h, i)  and  d(h, i)  will 

be  used  as  de- 

fined  in  Section  2. 

The  constraints  for  artery  ij  are  as  follows. 

For  the  artery 

as  a  whole: 

<  uiJ  <  (l/el^)z  , 

(LP4.1a) 

<  ulJ  <  d/eij)*  i 

(LP4.1b) 

for  each  signal  S^ 

on  the  artery: 

i 

j 

IJ  + 

wk  * 

‘1J  s  i  -  <k1J  . 

(LP4.2a) 

i 

i 

wlj  + 
k  + 

blJ  <  1  -  r  lJ  * 

—  k  * 

(LP4.2b) 
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for  each  pair  of  adjacent  siganls  on  the  artery,  say,  followed 

by  S-  outbound: 

X 

(wkij+  w'kij)  -  (wxiJ+  w/j)  +  d(k,5)  (uiJ+  Lij) 


(LP4.3) 


=  m(k,i)  -  (rk-rj)> 


m(k,jO  =  integer 


(LP4.4) 


4.2  Objective  Function.  In  a  progression  the  bandwidth  de¬ 
termines  the  maximum  length  (in  time)  of  a  platoon  that  can  pass 
unimpeded  down  the  street  when  the  prescribed  speeds  are  main¬ 
tained.  Therefore,  when  a  progression  system  is  appropriate,  the 
streets  (or  directions  on  streets)  should  presumably  be  allocated 
bandwidth  in  some  manner  according  to  the  flows  on  them.  Two 
devices  are  available  for  allocating  bandwidth  among  streets  (and 
directions):  the  objective  function  and’ the  constraints. 

Ideally  we  might  constrain  each  bandwidth  to  be  greater  than 
some  specified  size  needed  to  pass  a  known  flow.  Unfortunately, 
there  is  no  guarantee  that  the  result  will  be  feasible.  As  an 
alternative,  we  can  choose  some  important  bandwidth,  call  it  b°, 
and  maximize  this,  but  require  that  each  other  bandwidth  be 
greater  than  or  equal  to  a  specified  fraction  of  b°.  The  frac¬ 
tions  would  be  chosen  according  to  relative  flows. 

Such  constraints  are  not  sufficient  to  get  the  most  out  of 
the  system.  Sometimes  the  bandwidth  6n  an  artery  will  not  be  in 
conflict  with  bandwidth  on  other  arteries.  In  order  to  guarantee 
that  such  a  bandwidth  is  maximised,  every  bandwidth  must  be  in¬ 
cluded  in  the  objective  function.  The  less  important  bandwldths 

can  be  weighted  by  small  coefficients  and  b°  by  a  Urge  coefficient. 
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O  1  q 

To  formalize  these  ideas,  let  b  ,  b  ,...,b  be  the  band- 
wiJths  of  the  problem  (reindexed  here  for  convenience)  and  let 

a*  =  weight  assigned  to  the  i^  bandwidth;  i=0,...,q. 
k*  s  fraction  of  b°  to  be  guaranteed  to  the  i^  band¬ 
width;  i=l,...,  q. 

We  take  as  the  objective  of  the  network  program 
o  o  q  q 

max  a  b  +  . . .  +  a  b  ,  (LP4.5) 

as  constrained  by 

b1  >  k*  b°  ,  i=l, . . ,,q  (LP4.6) 

i 

Some  care  must  be  exercised  in  the  choice  of  the  k  so  as  not  to 

o 

let  an  unimportant  artery  limit  b  . 

4.3  Cycle  Constraints.  A  new  constraint  with  a  new  integer 
variable  must  be  introduced  whenever  several  arteries  intersect 
to  form  a  closed  loop  or,  as  we  shall  call  it,  a  cycle.  An  as 
example  of  a  cycle,  S^,  and  appear  at  the  arterial  in¬ 
tersections  of  a  cycle  in  Figure  6.  The  basic  reason  for  the  new 
constraints  lies  in  a  physical  requirement  of  simultaneity. 
Suppose  we  set  a  master  clock  to  aero  at  a  center  of  red  of  S^. 

If  we  proceed  down  artery  13  to  S^,  the  time  of  the  center  of  red 
there  is  fixed  by  the  progression  on  13.  As  a  result,  the  time 
of  the  center  of  red  for  along  the  direction  of  artery  35  is 
also  fixed  on  the  master  clock.  The  progression  on  35  then 
fixes  the  time  of  center  of  red  at  Sj.  Continuing  around  the 
cycle,  we  eventually  return  to  determine  the  time  of  a  center  of 
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red  at  Sp  This  time  must  be  an  integral  number  of  periods 
after  time  zero.  The  algebraic  statement  of  this  requirement 
constitutes  the  new  constraint. 

The  terminology  of  graph  theory  may  be  used  to  describe 
networks.  Signals  appear  at  nodes.  Lengths  of  street  between 
adjacent  signals  become  arcs  or  edges ♦  An  arc  has  a  direction 
and  is  usually  indicated  in  diagrams  by  an  arrow.  For  our 
purposes  the  direction  will  be  that  chosen  as  outbound  on  the 
street  segment.  An  edge  represents  the  street  segment  without 
regard  to  direction.  A  sequence  of  edges  that  forms  a  closed 
loop  is  called  a  cycle.  The  node  for  will  be  denoted  simply 
by  k.  The  arc  from  i  to  j  will  be  denoted  (i,j).  Edges  will 
be  referred  to  by  mentioning  the  corresponding  arcs  but  ignoring 
their  direction. 

Both  the  edge  and  arc  concepts  enter  the  formulation:  cycles 
are  formed  from  edges  but,  in  writing  down  equations,  it  is 
important  to  pick  a  direction  for  each  street  segment  and  use  the 
direction  consistently. 

The  formulation  of  the  constraint  for  a  single  cycle  pro¬ 
ceeds  as  follows:  Suppose  we  trace  around  a  cycle  and  find 
pieces  of  arteries  i^>  ^2^2,***^‘p^p  *****  oc<*el:,  Suppose 

that  the  signals  at  which  arteries  intersect  are  represented  by 
the  nodes,  k^  ^...kp  with  k^  at  the  intersection  of  arteries 
ipjp  and  i^,  k^  at  the  intersection  of  i^J^  and  k^J^,  etc.  The 
cycle  will  then  be  denoted  C(k^, . . . ,k^) . 
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Let  t  -  0  be  the  time  (measured  in  cycles)  of  a  center  of 


red  at  k^  along  artery  j ^ .  Then  t  ^(k^  k2>  is  the  time  of 
a  center  of  red  at  k^  along  artery  i  j  .  The  same  instant  is 
also  the  center  of  green  at  k2  along  artery  ^ 2^ 2 '  Therefore  a 
center  of  red  at  k^  along  occurs  at 

t  -  0<ki;  kp  +  1/2.* 


The  use  of  1/2  assumes  that  red  time  on  one  street  of  the  inter¬ 
section  coincides  with  green  time  on  the  cross  street  and  vice 
versa.  More  complicated  arrangements  are  possible  and  might 
require  replacing  the  1/2  by  some  other  fraction. 


Starting  from  this  time  we  can  add  0(k2^k^)  to  find  a  center  of 
red  at  along  artery  i^j^  •  Adding  1/2  we  turn  the  corner  and 
get  on  artery  i^j^  •  Proceeding  around  the  cycle,  we  eventually 
determine  the  time  at  which  a  center  of  red  occurs  at  k^  along 
artery  i^j^  .  Let 

c(k1,...,kp)  =  integer  variable  for  the  cycle  c^, . . . ,k^) . 
Then  the  cycle  constraint  is: 


c(k,,...,k  )  =  0(k,,kj+0(k  ,k,)+...+0(k  ,k  )+0(k  ,k)+(p/2).  (LP4.7) 

1  p  12  Z  J  p-1  p  pi 


Appropriate  expressions  for  0(i,j)  in  terms  of  other  variables  of 
the  program  are  given  by  (2.2)  or  (2.13),  possibly  with  some 
help  from  (2.4)  and  (2.5). 

Multiple  cycles  introduce  multiple  constraints.  The  number 
of  possible  cycles  is  greater  than  the  number  of  required  con¬ 
straints.  The  minimum  number  of:  constraints  can  be  developed  by 
tracing  over  the  network  as  follows;  Pick  a  starting  node.  Trace 
over  the  arcs  (irrespective  of  their  direction)  that  make  up  each 
artery  passing  through  the  starting  node.  Then,  for  each  node  reached 
by  the  tracing  so  far,  trace  over  the  arcs  on  the  arteries  that  pass 
through  the  node.  Now,  whenever  the  tracing  procedure  closes  a 
cycle,  form  a  constraint  from  the  cycle.  Continue  the  process 
until  no  further  arcs  can  bo  traced.  If  the  whole  network  has 
been  traced,  the  job  is  finished.  Otherwise,  the  network  breaks 
down  into  two  (or  more)  disconnected  sub-networks,  one  of  which 
has  Just  been  isolated.  Pick  another  starting  node  and  continue. 

The  tracing  process  works  because  it  identifies  which  signals 
bear  rigid  tlsw  relationships  to  the  starting  signal  by  virtue  of 
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being  connected  to  it  through  a  sequence  of  arterial  progressions. 
Whenever  a  cycle  is  closed  the  rigid  relationships  must  be  made 
compatible  by  a  constraint.. 

If  the  network  has  a  signal  at  every  intersection  of  arteries 
(i.e.  there  are  no  overpasses),  the  graph  of  the  system  is  planar. 
Then  the  number  of  cycle  constraints  equals  the  number  of  distinct 
areas  enclosed  by  arteries.  Figure  6  shows  a  planar  graph  re¬ 
quiring  two  cycle  constraints.  Cycle  constraints  are  not  neces¬ 
sarily  unique.  If  Figure  6,  for  example,  we  could  use  the  pair 
from  C(l,3,4,7)  and  C(7,4,5,6)  or  equally,  well  the  pair  from 
C(l,3,4,7)  and  C(l,3,5,6).  Algebraically,  either  of  these  pairs 
can  be  obtained  from  the  other. 

4.4  Red  Time  Constraints.  A  signal  at  the  intersection  of 
two  arteries  may  be  a  critical  signal  for  one  but  not  the  other. 
Then  a  shift  of  green  time  from  a  street  to  the  cross  street  may 
increase  bandwidth  on  one  artery  without  reducing  it  on  the  other. 
Thus,  the  red  time  at  arterial  intersections  may  sometimes  be  a 
useful  decision  variable.  The  variable  will  presumably  be  con¬ 
strained,  The  fraction  of  time  red  on  each  street  may  be  lisiited 
because  of  flow  requirements.  The  absolute  time  red  may  be  limited 
on  the  low  side  to  permit  pedestrian  crossing,  and  on  the  high  side 
to  placate  drivers. 

For  signals  at  which  red  time  is  a  decision  variable,  let 
B^,  «  lower  and  upper  limits  on  r^.  (cycles) 

Gk^,  ■  lower  and  upper  limits  on  r^T.  (seconds) 
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Then  the  variable  r^  is  constrained  by 

Ekij 

ij  ij 

-  rk  <  Fk 

(LP4.8a) 

Gkljz 

ij  ij 

<  rk  S  Hk  2 

(LP4.8b) 

4.5  Steps  in  Formulating  a  Network  Problem.  The  general 
mixed-integer  linear  program  for  the  network  problem  is  notation- 
ally  soraevliat  cumbersome  and  rather  than  write  the  program  out  we 
list  the  steps  that  generate  it. 

1.  Determine  which  streets  of  the  network  are  to  be 

arteries. 

2.  Set  up  the  objective  function  and  related  con¬ 
straints,  (LP4.5)  and  (LP4.6) 

3.  Set  up  the  period  constraints: 

(1/T2>  <  *  <  (1/T^  (LP4.9) 

4.  For  each  artery,  set  up  the  constraints, 

(LP4.1)  -  (LP4.4),  or,  if  preferred,  the  constraints  of  LP2  or 
LP3. 

5.  Set  up  cycle  constraints  using  (LP4.7)  and  the 
procedure  of  Section  4.3. 

6.  Add  red  time  constraints  (LP4.8)  for  any  red  times 
being  made  decision  variables. 

Network  problems  can  be  solved  by  tha  branch  and  bound  method?  of 
Section  3  with  a  few  modifications.  As  there,  the  integer 
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variables  are  introduced  one  at  a  time.  At  any  stage,  specific 
values  are  given  to  some  of  the  integer  variables  and  serve  to 
define  a  subset  of  solutions  to  the  network  problem.  An  upper 
bound  on  the  objective  function  for  solutions  in  the  subset  is 
found  by  solving  an  ordinary  linear  program  using  the  specific 
values  of  the  integer  variables. 

However,  in  the  network  problem,  the  full  n-signal  objective 
function  must  be  used  even  for  an  r-signal  problem.  Otherwise  the 
upper  bound  obtained  is  inappropriate.  As  a  result,  however, 
every  bandwidth  that  appears  in. .the  objective  function  must  have 
a  constraining  equation  (LP4.2)  in  each  linear  program  solved,  for 
otherwise  the  objective  function  will  have  an  unbounded  maximum. 
Therefore,  the  "r-signal"  problem  will  usually  include  constraints 
from  more  than  r  signals.  It  would  in  fact  be  permissible  to  use 
all  the  constraints  of  the  original  problem  (except  those  involving 
integer  variables  not  yet  assigned  specific  values)  in  every 
ordinary  linear  program.  With  the  exception  of  the  unboundedness 
problem  just  noted,  however,  the  constraints  associated  with  a 
signal  would  not  be  expected  to  affect  the  objective  function  very 
much  until  the  constraint  (LP4. 3)  containing  the  integer  variable 
is  added.  Thus,  adding  the  constraints  signal  by  signal  holds 
down  the  siae  of  the  linear  program  and  may  be  expected  to  save 
computation  time. 

The  order  of  adding  signals,  howaver,  will  affect  the  sise  of 
the  tree.  By  and  large,  troublesome  constraints  should  be  intro¬ 
duced  early.  Examples  would  be  constraints  for  signals  close 
together  and  constraints  for  cycles  that  have  ssmII  perimeters. 


4,6  Example.  Figure  6  displays  a  7  signal,  5  artery  prob¬ 
lem  which  has  been  solved  in  the  symmetric  case.  The  arteries 
are  13,  35,  47,  56,  and  16.  The  objective  function  used  is 


13  /l35  .16  56  47, 

b  +  .01  (b  +  b  +  b  +  b  ) , 


along  with  the  constraints 


b35>.5b‘3,  b16>  ,5b13,  b56>.5b13,  b4?  >  ,5b13 


Limits  on  period  are  50  and  100  seconds.  Speed  on  each  artery  is 

limited  by  14  and  16  meters/second  (31.3  and  35.8  mph) .  All  red 

13 

times  are  taken  to  be  .5  cycles  except  r^  ,  which  is  made  a  de- 

13 

cision  variable  subject  to:  .4  <  r^  <  .6  cycles  and 

13 

25  <  r^  T  <  50  seconds.  Distances  between  signals  are  shown  in 
Figure  6. 

The  complete  program  is  given  in  the  appendix  along  with  the 
solution  tree.  There  are  44  constraints  and  10  integer  variables. 
However,  in  the  branch  and  bound  solution  we  deal  with  42  con¬ 
straints  and  8  integer  variables.  The  reason  is  that  each  of  the 
two  cycle  constraints  is  introduced  by  limiting  the  values  that 
some  regular  Integer  variable  is  permitted  to  range  over. 

The  network  example  is  not  hard  to  solve.  15  linear  programs 

13 

were  required.  In  the  optimal  solution  the  bandwldths  are  b  =  .35, 
b35  ffl  .286,  b^  ■  .5,  b^  ■  .5,  and  b1^  >  .286  cycles.  The  period 
is  62.5  seconds.  The  speeds  are  v13  =  16,  v^5  ■  14,  v*7  •  16, 
v56  x  16,  and  ■  14  meters/second.  The  value  of  r^  *  .5 

cycles. 
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Figure  6.  A  network  of  7  signals.  In  this  network,  five  arteries;  13, 
33,  56,  16,  and  47,  have  keen  selected  for  consideration.  Two  cycle  con¬ 
straints  are  required,  arrows  Indicate  direction  chosen  as  outbound. 
Numbers  on  street  sagenU  are  distances  in  maters.  The  red-green  split 
has  been  made  a  dadaism  variable  at  7. 
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5.  Summary 

A  variety  of  maximal  bandwidth  problems  have  been  formulated 
as  mixed-integer  programs.  These  include  single  arterial  problems, 
more  complicated  arterial  problems  using  speed  and  period  as  de¬ 
cision  variables,  and  problems  involving  street  networks.  Some  of 
these  cases,  including  that  of  the  street  network,  do  not  appear  to 
have  been  put  in  any  kind  of  optimizing  format  before. 

Branch  and  bound  methods  are  given  for  solving  the  programs 
generated  and  examples  are  worked  out.  Although  it  is  not  known 
at  this  point  how  large  a  program  can  reasonably  be  solved  by  these 
methods,  it  seems  likely  that  problems  large  enough  to  be  of 
practical  interest  can  be  solved. 

★  *  *  ★  * 
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Appendix  1 


The  mixed' :  "  t  ^,er  linear  program  for  the  network  of  Figure  6 
is  given  be’ow  uv  t  e  symmetric  case. 

Find  blj,  z,  uiJ,  ,  m(k,*),  c  (i,j,M),  and  rfciJ  to 

L13  .  ,.35  .16  .56  47, 

max  b  +  .01(b  +  b  +  b  +  b  ) 

Subject  to: 


bandwidth  constraints: 

b35 

>  5b13 

(A.l) 

b16 

>  ,5b13 

(A.  2) 

b56 

>  -SO” 

(A.  3) 

b47 

>  -5b13 

(A.  4) 

period  constraints: 

.01 

VI 

N 

VI 

02 

(A.  5, 6) 

artery  13: 

,0625a  <  u13 

<  .0714a 

(A. 7,8) 

13 

V 

1 

+  b13  < 

.5 

(A.9) 

13 

W2 

+  b13  < 

.6 

(A. 10) 

13 

W1 

•  w213  +  200U13  =  . 5m(l,2)  -  .05 

(A. 11) 

• 

13 

W3 

+  b13  < 

.5 

(A. 12) 

*213-  w313  +  300a  3  .  . 5n(2, 3)  +.  05  (A. 13) 


i 


artery  35: 

•0625z  <  u  <  .07142 

(A. 14, 15) 

»335  +  b35  <  .5 

(A. 16) 

w“  +  b35  <  .5 

4  — 

(A. 17) 

+  150u  =  ,5m(3,4) 

(A.  18) 

35  .35 

Wj  +  b  <  .5 

(A. 19) 

35  35  35  c  ..  ev 

+  250u  =  .501(4,5) 

(A. 20) 

artery  56: 

.0625z  <  u56  <  .0714z 

(A.  21, 22) 

56  56 

+  b  <  .5 

(A.  23) 

56  l56  ^  _ 

wg  +  b  <  .5 

(A. 24) 

56  56  56  _ 

"  wg  +  500u  =  .5ra(5,6) 

(A. 25) 

artery  47: 

47 

.0625*  <  u  <  .0714* 

(A. 26, 27) 

w4  +  b  <  .5 

(A. 28) 

Wj  +  b  <  l  -  ty 

(A. 29) 

Red  tine 

47  ^  16  , 

r  +  r?  =1 

(A.  30) 

constraints 

.  16 
•4  £  r7  -  *6 

(A.  31, 32) 

16 

25s  <  r?  <  50* 

(A.  33, 34) 

47  47  47  47 

w4  -  »?  +  500u  .  .5n(4,7)-.5(.5-r7  ) 

(A.  35) 

39- 


artery  16; 


,0625z  <  u  <  ,0714z 

16  16 

w1  +  b  <  -5 

16  .16  ,  16 

Vy  +  b  <  l  -  r^ 


(A.  36, 37) 
(A. 38) 


(A. 39) 


Wl16-  w?16+  150u16=  .5ra(l,7)T.5(.5-r716)  (A. 40) 


v“  +  b16  <  .5 
0 


(A. 41) 


w_16-  w,16+  250u16=  .5m(7,6)-.5(r,l6-.5)  (A. 42) 
/  o  7 


cycle  C(l,3,4,7): 


m(l,2)  +  o(2,3)  +• 0(3,4), +  »(4,7)  -  «(1,7) 


2c(l, 3,4,7) 


(A.  43) 


Cycle  C(4,5,6,7) ; 


o(4,5)  +  o(5,6)  -  o(7,6)  -  o(4,7)  = 
2c(4,5,6,7) 


(A.  44) 


«(M)i  c(i,j,k,je)  integers 
wiJ,  biJ  >  0 


In  using  the  branch  end  bound  algoritho,  the  constraints  were 
introduced  in  the  following  blocks,  each  terminating  with  an  integer 
variable  constraint:*  (1)  A.1-A110,,  A. 16,  A. 23,  A. 28,  A. 38,  A. 11  . 
(2)  A.  12,  A.  13  .  (3)  A.  14,  A.  15,  A.  17,  A.  18  .  (4)  A. 26,  A. 27, 
A.29-A.35  .  (5)  A. 36,  A. 37,  A. 39,  A.40,  A. 43. .  (6)  A. 19,  A.20  . 
(7)  A.  41,  A. 42  .  (8)  A. 21,  A. 22,  A. 24,  A. 25,  A.44  . 


m(i,2) 


m(2,3) 


m(3,4) 


m(4,7) 


nul,?>  ' 


m(4,5) 


m(7,6) 


m(5,6) 


Figure  7.  Tree  for  network  example.  Labels  on  the  nodes  are  the  calculated 
upper  bounds  on  the  objective  function.  The  objective  function  is  a  weighted 
sum  of  arterial  bandwidths.  Numbers  inside  the  nodes  are  values  for  the  m. 
shown  at  right.  Cycle  constraints  have  been  used  to  limit  the  possible 
choices  of  m(l,7)  and  m(5,6). 
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Appendix  2 


The  mixed-integer  program  for  the  10  signal  problem  is  given 
below. 


Find  b,  z,  t^  m^  to 

max  b 


Subject  to: 

.0182  <  z  <  .0133  (A2.1,2) . 

wx  +  b  <  .53  (A2.3) 

9.4z  <  t.  <  12. 5z  (A2.4,5) 

+  b  <  .60  (A2 . 6) 

*i  -  w2  +  t  =  (l/2)o»l  -  .035  (A2.7) 

11. 9z  <  t2  <  15.9*  (A2.8,9) 

-2.03z  <  .787t2  -  t  <  2.03z  (A2.10,il) 

w3  +  b  <  .60  (A2.12) 

W2  "  W3  +  C2  *  (l/2)n>2  (A2.13) 

18.7*  <  t3  <  25.0*  (A2.14, 15) 

-2.58*  <  .636t3  -  t2  <  2.58*  (A2.16,17) 

w4  +  b  <  .53  (A2.18) 


w3  ‘  W4  +  C3 


<l/2)m3  +  .035 


(A2.19) 


f 


•7 


11.9*  <  t.  <  15. 9z 
—  4  — 

-4.05z  <  1 . 571 t .  -  t,  <  4.05z 
-  4  3- 

+  b  <  .52 

w.  "  w  +  t  =  (l/2)m  +  .005 

4  5  4  4 


(A2.20,21) 
(A2. 22, 23) 
(A2.24) 
(A2.25) 


13. 6z  <  t5  <  18. 2z  (A2. 26, 27) 

-2.58z  <  .873t5  -  t4  <  2.58z  (A2. 28,29) 

wfi  +  b  <  .58  (A2„30) 

wc  -  w,  +  t  =  (l/2)m  -  .03  (A2.31) 

5  o  5  5 

11. lz  <  <  14. 8z  (A2. 32, 33) 

-2.95z  <  1.231t6  -  t$  <  2.95z  (A2.34,35) 

w?  +  b  <  .60  (A2.36) 

v6  -  v?  +  t  =  (l/2)m6  -  .01  (A2.37) 

6.8z  <  t?  <  9,  lz  (A2. 38, 39) 

-2.40z  <  1.624t?  -  t6  <  2.40z  (A2.40,4i) 

wg  +  b  <  .60  (A2.42) 

w?  -  wg  +  t?  =  (l/2)m?  (A2.43) 


*»3' 


11. 9z  <  tg  <  15. 9z 


(A2.44,45) 


-1.48z  <  .573tg  -  t?  <  1.48z 

+  b  <  .60 

w8  "  W9  +  fc8  =  (1/2)m8 


(A2.46,4? 

(A2.48) 

(A2.49) 


7. 7z  <  t9  '<  10. 2z 


-2.58z  <  1.555tq  -  t  <  2.58z 


w10  +  b  <  -58 


w9  ’  W10  +  C9  =  <1/2>m9  +  *01 


wi>  b  >  0 


(A2.50,51) 
(A2. 52,53) 
(A2.S4) 
(A2.55) 


The  above  data  is  arranged  in  blocks.  The  first  (r-1) 
blocks  represent  an  r-signal  problem  for  signals  sj>“,;sr* 
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